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ABSTRACT. 

In the Sherrington-Kirkpatrick mean field model for spin glasses, we show that the quen- 
ched average of the free energy can be expressed through a couple of functional order 
parameters, in a form very similar to the one found in the frame of the replica symme- 
try breaking method. The functional order parameters are implicitely given in terms of 
fluctuations of thermodynamic variables. 

Under the assumption that the two order parameters can be chosen to be the same, in 
the thermodynamic limit, it is shown that the Parisi free energy is a rigorous upper bound 
for the free energy of the model. 
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Let us introduce the partition function Z^{fi, J) and the free energy Fn(/3, J) for the 
Sherrington-Kirkpatrick mean field spin glass model [1,2] in the form 



Z N ((3,J) = exp( JL_ J2 JijWj) = exp(-PF N (P, J)). (1) 

CT1...CTJV 

The cr's are Ising spins describing a generic configuration of the system 

a: {1,2,..., N} 3 i-KneZ 2 = {-1,1}. (2) 

For each of the N(N — l)/2 couples of sites i ^ j, over which the sum j) 

runs, we have introduced independent random variables = Jji, i ^ j, identically 
distributed , called quenched variables. The cr's are mesoscopic random variables subject 
to thermodynamic equilibrium. The J's do not participate to thermodynamic equilibrium, 
but act as a kind of random environment on the cr's. For the sake of simplicity, we assume 
that the J's have unit Gaussian distribution with 

E(J tJ )=0, £(4) = 1, (3) 

where E denotes averages with respect to the J variables. The parameter (5 is the inverse 
temperature in proper units. 

We are interested in the expression of the thermodynamic limit N — > oo for the free 
energy per spin, averaged over the external noise (quenched average), 

lim W-^logZjvCftJ)). (4) 

Let us introduce the marginal free energy, i.e. the increment in the free energy when an 
additional (iV + l)th spin is added to a system of iV spins, at the same inverse temperature 

-f3(F N+1 (/3, J) - F N (p, J)) = \ogZ N+1 (p, J) - logZ N (p, J). (5) 

Then we have 

Proposition 1. The quenched average of the marginal free energy and the free energy 
per spin can be expressed in the following form 

E(\ogZ N+1 (p, J)) - E(\ogZ N (p, J)) = lo 9 2 + <M/3) - <f> N (p), (6) 
(N + l)- 1 E(logZ N+1 ((3,J)) = log2 + ^ N (p)-0 N ((3), (7) 
ip N (/3) = E log u N (cosh-=y2 J^), (8) 

(j> N (/3) =E\ogu N (exp ^ N ^ N _ ^ Yl Ji3 a i a j)i ( 9 ) 
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iP N (/3) = (N + 1)- 1 J2 <Pn(P) = (N + 1)- 1 J2 4>k(P). (10) 



K=0 K=0 

Here a; at is the Boltzmann state with Boltzmannfaktor as in (1) with y/N — 1 replaced 
by \/N, and the Jj's and J^'s, i, j = 1, 2, . . . , N, i ^ j, are N and N(N — l)/2, respectively, 
independent random variables, with unit Gaussian distribution. We call Jij the stale noise, 
and Ji and Jij the fresh noise. 

For the proof we can write 

1...N 1...N 



E (log Z N+l (f3, J)) = Elog V* exp(^= V J ij a l a j + V* -/iiv+i 0-^+1) = 

(T1...(TJV + 1 (ijj) t 

= log2 + £ log wjv( cosh -^=J^ J i£7i ) + £log exp(-^=^ JijWj), (H) 



/A 

(T1...(TJV + 1 

> exp(^ 

(71...CTJV 

where we have explicitely performed the sum over o~n+i and have called Ji the old JiN+i- 
Let us now consider 

E\ogZ N ((3,J) = E\og exp( P— YlJiiWj)' ( 12 ) 

(71...CTJV 

By introducing a fresh set of independent noise Jy, with the same normalization as 
in (3), we can replace Ji j / \J N — 1 with the stochastically equivalent sum J^j/y/N + 
Jij/ a/ A 7 " (A 7 " — 1), in fact the two random variables have the same mean and tha same 
covariance. Therefore, we have 

E\ogZ N (/3, J) = Elog J2 exp ^7^ Yl JijWj + /N L _ /Z Ji3°i°3) = 



/; log (exp(— ===== J2 JijWj)) +Elog e ^77^ Yl ^i^i)' ( 13 ) 



and (6) follows. Now we can write (6) for a generic K, and sum from K = to K = N . 
With the obvious notations 

£logZoC6U) = 0, Vo(/3) = 0, 0o(^) = O, (14) 

we immediately have (7). 

Useful information on the functions ^jv(^), ip N (/3), 4>n(P), 4>n(P) * s given by 
Theorem 2. The following bounds hold 

1 logcosh(/3z) < <M/3), ^(/3) < /3 2 /2, (15) 
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0<<M/3), <f> N (P) </3 2 /4, 



(16) 



where d/j,(z) = exp(—z 2 /2) dz/ y/2~7v is the unit Gaussian distribution. 

For iPn(/3) the proof has been given in [3]. It is based on either annealing the E 
averages (upper bound), or quenching the ui averages (lower bound). The bound for 
4>n{P) follows easily from the definition (10). In the same way one proves (16). 

Let us now introduce the convex set X of functional order parameters of the type 

x , [0,1] 9^x^)6 [0,1], (17) 

with the L 1 (<iq i ) distance norm. We induce on X a partial ordering, by defining x < x if 
x(q) < x(q), for all < q < 1, and introduce the extremal order parameters xo(q) = and 
Xi(q) = 1, such that for any x we have Xo(q) < x(q) < Xi(q). 

For each x in X, and (3 > 0, let us define the function with values f(q,y;x, (3), 
< q < 1, y E R, as the solution of the nonlinear antiparabolic equation 

d q f+l(f" + x(q)f' 2 ) =0, (18) 

with final condition 

/(l,y;x,0) = logcosh(/fy). (19) 

In (18), /' = d y f and /" = d 2 y f. 

As a shorthand notation, for each x in X, and /3 > 0, we define at q = 0, y = 

f(x,(l) = f(0,0;x,(3). (20) 

In Ref. [3], we have shown that (18,19) arise in a very natural way as a result of exact 
corrections to the annealing approximation logi?a;(. . .) in the evaluation of the quenched 
average in (8). 

The following theorem summarizes some important properties [3] of f(x,/3). 

Theorem 3. The function f(x,/3) is monotone in x, i.e. x > x implies f(x,/3) > f(x,/3). 
Moreover, the following bounds hold 

f(x ,P) = J logcosh(^) dfJL(z) < f(x,f3) < (5 2 /2 = f(x 1 ,f3). (21) 

In [3], we have also proven the following representation theorem. 

Theorem 4. There exists a nonempty hypersurface T<n{/3) in X, such that, for any x G X 
and / solution of (18,19), we have the following representation 

iP N (P) = f(x,(3). (22) 

Any family of functional order parameters, x e , depending continuously in the L 1 norm 
on the variable e, < e < 1, with x$ = 0, and xi = 1, and nondecreasing in e must 
necessarily cross Ejv(/3) for some value of the variable e (we say that Ejv(/3) has the 
monotone intersection property). 
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By using the same method, we can easily prove the following easy generalization. 
Theorem 5. The average ip N (/3), defined in (10), admits also a representation 

^v(/3) =/(x,/3), (23) 

for x on some hypersurface Eat(/?) (which will be in general sligthly different from the 
hypersurface appearing in the previous (22)). 

The method of Ref. [3] allows to give implicit expressions for the elements of Ejv(/3) 
in terms of fluctuations, but the very existence of Ejv(/3), with the monotone intersection 
property, follows from a very simple argument. In fact, from the bounds (15) and (21), and 
the monotonicity of f(x,/3) in x, given by Theorem 3, we immediately have the existence 
of a nonempty £jv(/3). 

Similar representation formulae hold for 4>n{P) and 4> N (f3). 

Theorem 6. There exist nonempty convex linear sets Yi' N ((3) and Ejv(/?) m X, such that 

1 r 1 

cj> N (P) or cj> N (p) = -P 2 J q x(q) dq, (24) 

for any x E T,' N (/3), or x, E Sjv(^), respectively. 

The proof follows from a simple cumulant expression. Let us introduce the interpo- 
lating parameter q, < q < 1, and define 

(j){q) = Elogu N (exp ^== _ ^ J2 ^i a ^i)^ ( 25 ) 



so that (p(0) = and 0(1) = 4>n{P), as defined in (9). Let us take the derivative 

<J ? >[q) = ~7WtW=t{ 

Then we can exploit the general integration by parts formula 



dq <t>(<l) = y / N ( N _ 1 ) J2 E (^ UJ N 1 ( e M- ■ -^N^jexpi- ••)))• ( 26 ) 



and obtain 



where 



E(J ij F(J)) = E(^F(J)), (27) 

OJij 



±J>(q) = (28) 



_ yNiwi ex P T^rry Jaw) 

x(q) = l >E( ^— )• 29 
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Clearly, we have the inequality 

< x(q) < 1. 



(30) 



By integrating (28) on dq we have the representation (24) for 4>n{(5). This shows that 
T,' N ((3) is nonempty, because x is explicitely defined by (29). Of course, all functional order 
parameters, which give the same value for the integral in (24), are acceptable. This is how 
the convex linear set T,' N ((3) arises. Also in this case we have the monotone intersection 
property. The representation (24) for </> N (/3) follows easily from the definition (10). 

Let us also explicitely remark that the representations given in Theorems 4,5,6 hold 
for any even state ojn, not necessarily as that arising in (8,9). Of course, the involved 
hypersurfaces do depend on the particular u>n- 

By collecting all results of Theorems 5 and 6, and the definition (7), we have the 
following basic representation theorem for the quenched average of the free energy per 
spin 

Theorem 7. There exist nonempty hypersurfaces T, N (/3) and £jv(/3) in X, such that 

1 r 1 

(N+l)- 1 E{\ogZ N+1 (p,J))=\og2 + f(x,p)--p 2 j^ qx(q) dq, (31) 

for any x G £at(/3) and x, G Ejv(/3). Elements of these two hypersurfaces can be expressed 
implicitely in terms of fluctuations. 

The representation (31) is equivalent and complementary to the representation given 
in [2], which involves the order parameter x for different values of (3. Here two order 
parameters are involved, but at the same value of (3. 

This representation is very similar to that found in the frame of the replica symmetry 
breaking method, with Parisi Ansatz [2], where the two order parameters are considered 
to be the same, at least in the thermodynamic limit. Therefore, we are led to explore the 
consequences of the following 

Assumption 8. Let a^^P) be the L 1 distance between the hypersurfaces £at(/3) and 
x G £ N (0) 

a N (0) = inf f \x{q) - x{q)\ dq, x G Z N {/3), x G £ N (0), (32) 
Jo 

and assume 

lim a N ((3) = 0. (33) 

AT— s-oo 

Let us also define the Parisi free energy fp(/3) at inverse temperature /3 as 

i r 1 

-pfp(/3) = inf (log 2 + f(x,p) --P 2 q x(q) dq). (34) 

Then we have 

Proposition 9. Under the stated assumption, in the thermodynamic limit, we have 

lim inf (N + l)~ 1 E(logZ N+ i((3, J)) > -f3f P (f3). (35) 
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The proof is immediate. In fact, for the r.h.s. of (31) we have 
log 2 + f(x, 0) - \p 2 J q x(q) dq + ^(3 2 J q (x(q) - %)) dq > -(3f P ((3) - a N (0), (36) 

and the result follows by taking the limit N — > oo. 

Therefore, the Parisi free energy, with these assumptions, is proven to be at least a 
rigorous upper bound for the infinite volume limit of the free energy of the model. 

In a forthcoming paper [4], we show that there is good evidence, not a definite math- 
ematical proof as yet, that the two order parameters in (31) can be taken the same, in the 
thermodynamic limit, and moreover that the Parisi free energy is the true free energy, and 
not only an upper bound. 
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